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A  multi-objective  parameter  identification  method  for  modeling  of  Li-ion  battery  performance  is  pre¬ 
sented.  Terminal  voltage  and  surface  temperature  curves  at  15  °C  and  30  °C  are  used  as  four  identification 
objectives.  The  Pareto  fronts  of  two  types  of  Li-ion  battery  are  obtained  using  the  modified  multi¬ 
objective  genetic  algorithm  NSGA-II  and  the  final  identification  results  are  selected  using  the  multiple 
criteria  decision  making  method  TOPSIS.  The  simulated  data  using  the  final  identification  results  are  in 
good  agreement  with  experimental  data  under  a  range  of  operating  conditions.  The  validation  results 
demonstrate  that  the  modified  NSGA-II  and  TOPSIS  algorithms  can  be  used  as  robust  and  reliable  tools 
for  identifying  parameters  of  multi-physics  models  for  many  types  of  Li-ion  batteries. 
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1.  Introduction 

In  recent  years,  Lithium-ion  batteries  have  been  widely  used  for 
electric  vehicles  and  energy  storage  due  to  their  relatively  high 
energy  and  power  density.  A  key  challenge  facing  battery  manu¬ 
facturers  is  the  limited  durability  of  the  devices  under  demanding 
operating  conditions.  As  a  result,  characterization  of  battery 
degradation  using  material  and  electrochemical  parameters  is  a 
popular  topic  and  has  attracted  many  researchers.  X.  Han  et  al.  [1] 
and  Q.  Zhang  et  al.  [2]  used  the  stoichiometric  number  of  the 
electrode  material  to  indicate  the  stages  of  capacity  fade.  A.P. 
Schmidt  et  al.  investigated  the  relationship  between  the  volume 
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fraction  of  active  material  in  electrode  (fs),  the  ionic  conductivity  of 
the  electrolyte  (/ce)  and  the  cycle  number  during  aging,  and  used 
them  as  characteristics  to  evaluate  the  state  of  health  (SOH)  of 
batteries  [3  .  V.  Ramadesigan  et  al.  found  that  both  the  electro¬ 
chemical  reaction  rate  (/<s,a)  And  the  Li+  diffusion  coefficient  (Ds,a)  in 
the  anode  degrade  with  cycle  number  following  a  power-law  [4  .  R. 
Fu  treated  some  parameters  of  an  electrochemical  thermal  model 
as  degradation  parameters  to  analyze  the  degradation  effects 
considering  side  reactions,  including  the  volume  fraction  of  active 
material  in  electrode,  the  solid  electrolyte  interface  (SEI)  resistance 
(fifilm)  and  the  diffusion  coefficient  of  electrolyte  (De)  [5].  These 
papers  identified  some  isolated  relationships  between  a  limited 
number  of  parameters  and  battery  health.  Our  research  pro¬ 
gramme,  by  contrast,  is  focused  on  the  investigation  of  a  much 
larger  set  of  model  parameters  in  different  stages  of  aging,  which 
requires  an  accurate  electrochemical  model  and  a  nondestructive 
method  for  parameter  identification. 

Several  electrochemical  models  have  been  presented  in  the 
published  literature.  M.  Doyle  et  al.  proposed  the  P2D  model  [6,7], 
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which  addressed  the  complex  mechanisms  of  physical  and  chem¬ 
ical  processes  such  as  diffusion,  transport  of  ions,  ohmic  phenom¬ 
ena  and  electrochemical  reaction,  using  a  group  of  partial 
differential  equations  and  some  algebraic  equations.  The  Single 
Particle  Model  (SPM)  has  a  simpler  form  and  fewer  parameters 
compared  to  the  P2D  model;  it  can  predict  diffusion  and  interca¬ 
lation  processes  but  the  local  current  distribution  is  neglected,  so  its 
use  is  only  valid  at  low  current  densities  [8].  The  Extend-SPM  model 
takes  into  account  the  effects  of  electrolyte  concentration  and 
electrochemical  reaction  distribution  [9  ,  but  its  predictive  capa¬ 
bility  is  still  inferior  to  the  P2D  model.  In  addition,  since  thermal 
behavior  is  not  included  in  the  Extend-SPM  model,  it  is  impossible 
to  identify  temperature  related  parameters.  The  electro-thermal 
model  is  a  P2D-based  model,  which  includes  some  thermal 
behavior  such  as  energy  conservation,  heat  generation  and  ex¬ 
change,  and  can  accurately  predict  a  range  of  phenomena  occurring 
in  the  battery  10].  P2D-based  models  have  a  large  set  of  parame¬ 
ters  with  corresponding  physical  meaning,  which  makes  them  very 
suitable  for  the  study  of  degradation  mechanisms.  Recently,  the 
experimentally  validated  simulation  results  of  P2D-based  models 
with  accurate  parameters  have  also  been  used  for  capacity  fading 
research,  such  as  the  degradation  behavior  caused  by  Li-plating  in 
low  temperature  charging  process  [11]  or  elementary  reaction- 
based  SEI  layer  generation  under  different  operating  conditions 
[12]. 

Because  of  the  high  complexity  and  nonlinearity  of  P2D-based 
models,  analytical  solutions  do  not  exist.  In  addition,  the  number 
of  parameters  is  large  (27  in  this  work)  so  parameter  identification 
is  very  challenging  using  conventional  approaches.  Genetic  algo¬ 
rithms  (GAs)  are  a  class  of  stochastic  algorithms  well-suited  for 
large-scale  optimization.  They  have  been  successfully  applied  to 
parameter  identification  for  equivalent  circuit  models  [13  and  P2D 
models  [14]  of  Lithium-ion  batteries.  The  rationale  of  the  GA  pro¬ 
cedure  is  to  simulate  the  process  of  natural  selection  and  mutation 
and  determine  the  optimal  solution  to  minimize  the  errors  be¬ 
tween  simulated  data  and  experimental  data. 

In  previous  studies,  the  constant-current  test  4]  or  the  dynamic 
charge/discharge  test  [14]  at  a  certain  ambient  temperature  was 
used  for  parameter  identification.  The  results  are  only  available  for 
certain  identification  conditions,  while  under  other  conditions 
(especially  at  other  ambient  temperatures)  the  identified  results 
may  not  predict  the  battery  performance  accurately.  The  multi¬ 
objective  optimization  approach  has  been  successfully  applied  to 
battery  modeling  by  previous  authors.  P.  Kumar  et  al.  15]  used 
Nondominated  sorting  Biologically  Motivated  Genetic  Algorithm 
(NBGA)  and  J.  Brand  et  al.  16,17]  used  Nondominated  Sorting  Ge¬ 
netic  Algorithm  (NSGA-II)  for  parameter  identification,  treating  the 
charge/discharge  data  at  different  currents  as  multi-objectives. 
Their  results  demonstrated  that  this  method  is  applicable  to 
many  types  of  batteries.  However,  the  method  has  yet  to  be  applied 
to  P2D-based  models,  and  the  surface  temperature  has  not  been 
included  as  an  identification  objective.  In  contrast  to  a  traditional 
genetic  algorithm,  the  multi-objective  genetic  algorithm  provides  a 
set  of  nondominated  solutions  instead  of  a  certain  optimal  solution 
as  an  identification  result.  There  remains  an  uncertainty  in  the 
literature  over  how  to  decide  the  final  identification  result  [15—17]. 

In  this  article,  the  parameters  of  our  multi-physics  model  are 
identified  using  a  multi-objective  genetic  algorithm;  the  remainder 
of  the  paper  is  organized  as  follows:  Section  2  outlines  the  multi¬ 
physics  model  of  a  Lithium-ion  battery.  Section  3  summarizes  the 
experimental  setup  and  test  data.  Section  4  presents  the  principle 
and  method  of  parameter  identification  using  the  NSGA-II  algo¬ 
rithm.  Section  5  provides  the  final  identification  result  and  presents 
validation  studies  by  comparing  the  simulated  curves  with  the 


experimental  data  under  a  range  of  operating  conditions.  Conclu¬ 
sions  of  this  article  are  presented  in  the  final  section. 

2.  Multi-physics  model 

The  multi-physics  model  consists  three  parts:  the  electro¬ 
chemical  model,  the  thermal  model  and  the  radial  heat  distribution 
model. 

2.1.  Electrochemical  model 

The  electrochemical  part  of  the  model  is  described  by  a  con¬ 
ventional  P2D  model.  Fig.  1  shows  a  schematic  of  the  P2D  model  for 
a  Lithium-ion  battery.  It  consists  of  two  current  collectors,  a 
negative  electrode  (anode),  a  separator  and  a  positive  electrode 
(cathode).  Both  electrodes  and  separator  are  porous  structures.  Two 
inner  boundaries  (anode/separator  interface  2  and  separator/ 
cathode  interface  3)  and  two  external  boundaries  (Cu/anode 
interface  1  and  cathode/Al  interface  4)  are  also  shown. 

The  physical  and  chemical  processes  such  as  ion  diffusion  and 
migration  and  electrochemical  kinetics  are  described  by  several 
partial  differential  equations  and  algebraic  equations  in  Table  1. 

In  Eq.  (6),  E0cv  is  the  equilibrium  potential,  which  is  a  function  of 
the  solid  phase  Li+  concentration  at  the  particle  surface.  The 
reference  E0 cv  of  the  graphite/LiFePCH  system  is  obtained  from 
Ref.  [18  ,  and  the  reference  E0 cv  of  the  MCMB/LiCo02  system  is 
obtained  from  Ref.  [8  .  The  meanings  of  the  symbols  are  explained 
in  detail  at  the  end  of  this  article. 

2.2.  Thermal  model 

Thermal  phenomena  such  as  energy  conservation,  heat  gener¬ 
ation  and  exchange  also  occur  during  the  charge/discharge  process, 
and  they  can  be  simulated  by  the  P2D  model,  as  presented  in 
Table  2.  The  heat  generation  occurs  in  the  spiral  roll,  including  the 
heat  of  electrochemical  reaction,  entropic  heat  and  ohmic  heating. 
Heat  exchange,  including  convection  and  radiation,  only  occurs  at 
the  surface  of  the  cylindrical  battery. 

Some  parameters  in  the  P2D  model,  such  as  solid  phase  diffu¬ 
sion  coefficients  Ds>a  and  DS)C,  solution  phase  conductivity  /ce,  solu¬ 
tion  phase  diffusion  coefficient  De  and  electrochemical  reaction  rate 
/<s,a  and  /<s,c,  are  coupled  with  the  battery  temperature.  These  pa¬ 
rameters  can  be  updated  by  the  Arrhenius  law.  In  Eq.  (15),  Xj  rep¬ 
resents  the  parameters  mentioned  above,  Xi>ref  is  the  parameter 


Fig.  1.  Schematic  of  P2D  model  (discharge  process). 
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Table  1 

Governing  equations  and  boundary  conditions  of  the  P2D  model. 
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Mechanisms 
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Boundary  condition 


Solid  phase  diffusion 

Electrolyte  phase  diffusion 
Solid  phase  charge  balance 

Electrolyte  phase  charge  balance 

Electrochemical  kinetics 
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value  at  reference  temperature  Tref,  Ei  denotes  the  corresponding 
activation  energy  of  parameter  Xh  respectively.  In  addition,  the 
equilibrium  potential  depends  on  battery  temperature  according  to 
the  Nernst  equation  (Eq.  (16)).  An  electrochemical-thermal 
coupled  model  is  presented  by  feeding  back  the  relationship  be¬ 
tween  temperature  dependent  parameters  and  thermal  effects  to 
the  P2D  model.  There  is  a  well  established  literature  on  the  P2D 
model  and  the  electrochemical-thermal  model,  the  reader  is 
referred  to  [6-8,18-20]  for  details. 

In  Eqs.  (11)  and  (16),  dE0cv/dT  is  the  entropy  coefficient,  and 
obtained  from  Ref.  [19  for  graphite/LiFePOq  system  and  Ref.  [8]  for 
MCMB/LiCo02  system. 


2.3.  Thermal  distribution  of  cylindrical  batteries 

A  typical  cylindrical  Li-ion  battery  is  a  spiral-rolled  cell  as  shown 
in  Fig.  2.  Thin  layers  of  anode,  cathode,  separator  and  current  col¬ 
lector  are  stacked  up  as  Section  B-B  and  rolled  up  on  central 
mandrel  as  Section  A-A,  and  then  inserted  into  a  can.  The  gaps  are 
filled  with  liquid  electrolyte. 

The  temperature  distribution  in  the  radial  direction  of  a  cylin¬ 
drical  battery  is  more  significant  than  that  in  the  axial  direction  due 
to  the  lower  thermal  conductivity  21  .  For  a  26650  type  cylindrical 
battery,  the  radius  is  0.013  m,  the  heat-transfer  coefficient  of  the 
lateral  surface  is  -10-50  W  m“2  K-1,  and  the  radial  thermal  con¬ 
ductivity  is  about  0.15  W  itT1  K_1,  so  the  Biot  number  is  greater 
than  0.1.  The  temperature  gradient  in  the  radial  direction  must 
therefore  be  taken  into  account,  otherwise  the  simulated  temper¬ 
ature  will  be  very  different  from  the  experimental  value  measured 
at  the  surface  of  the  battery  [22  .  The  heat  conduction  term  in  Eq. 
(9)  can  be  expressed  in  one-dimensional  by  Fourier's  Law: 


-V(ArVT)  =  Ar 
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The  temperature  gradient  at  the  center  of  the  battery  is  zero, 
while  at  the  battery  surface  it  is  given  by  the  heat  exchange  rate  at 
the  surface: 
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where  Rceii  is  the  radius  of  the  cylindrical  battery  and  the  heat 
exchange  rate  is  q  is  given  by  q  =  qc  +  qr,  where  qc  and  qv  are  the 
heat  exchange  rates  due  to  conduction  and  radiation  respec¬ 
tively.  The  temperature  gradient  is  continuous  at  the  interfaces 
of  spiral  roll,  gap  and  can.  The  temperature  of  the  can  is  defined 
as  the  shell  temperature  Ts h  which  can  be  measured  from 
experiment. 

A  thermal  impedance  model  has  been  developed  to  simulate  the 
thermal  distribution  in  the  radial  direction  of  cylindrical  Lithium- 
ion  batteries  by  solving  Eq.  (17),  as  detailed  in  our  previous  work 
[23].  In  contrast  to  other  studies  [14,18,19],  the  surface  temperature 
is  a  model  output,  facilitating  its  use  as  an  optimization  objective  in 
addition  to  the  terminal  voltage. 


3.  Experiments 

Two  types  of  Lithium-ion  battery  are  investigated  in  this  article; 
one  is  a  2.3  Ah  LiFePOq  battery  (A123  system  ANR26650mlA,  USA), 
and  the  other  is  a  1.8  Ah  LiCo02  battery  (Lishen  LS18650BC,  China). 
Before  testing,  the  batteries  were  charged  with  standard  constant 
current-constant  voltage  (CC-CV)  strategy  -  the  CC  rate  was  C/3 
and  the  CV  voltage  was  3.6  V  for  the  LiFePOq  battery  and  4.2  V  for 
the  LiCo02  battery,  and  the  cutoff  current  for  the  CV  process  was  C/ 
50.  After  a  30  min  relaxation,  simplified  dynamic  stress  tests  (DSTs) 
were  conducted  in  an  environmental  chamber  (Partner,  China)  at 
Tam  =  15  °C  and  Tam  =  30  °C,  respectively.  Experimental  Uapp  data 
were  acquired  by  a  Battery  Testing  System  (Neware  BTS-5V-6A, 
China)  every  10  s,  while  the  experimental  Ts h  data  were  acquired 
simultaneously  by  the  auxiliary  equipment  (Neware  BTS-AUX, 
China)  using  thermistors  as  sensors.  The  temperature  sensors 
were  attached  to  the  shell  of  the  batteries  using  thermally 
conductive  silicone.  The  simplified  DST  profiles  of  these  two  bat¬ 
teries  are  shown  in  Fig.  3. 
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Table  2 

Equations  of  the  thermal  model. 


Mechanisms  Equation 
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For  the  LiFeP04  battery,  there  are  two  discharging  rates  (1C  and 
2C)  and  one  charging  rate  (1C)  in  the  test,  while  for  UC0O2,  there 
are  three  discharging  rates  (C/4,  C/2  and  1C)  and  two  charging  rates 
(C/4  and  C/2).  There  is  a  trade-off  between  increasing  the  number  of 
charge/discharge  rates  to  make  the  identification  result  more  ac¬ 
curate  and  robust,  at  the  expense  of  test  complexity  and  compu¬ 
tation  time. 


The  area  covered  by  the  Pareto  front  is  called  the  hyper-volume. 
Any  point  inside  the  hyper-volume  is  considered  to  be  dominated 
by  at  least  one  solution  in  the  Pareto  front.  The  Pareto  front  itself  is 
nondominated,  means  there  is  no  solution  in  this  set  that  can 
optimally  satisfy  all  the  objectives  16]. 

In  the  standard  NSGA-II  process,  a  set  of  parameters  is  consid¬ 
ered  to  be  an  individual.  Firstly,  a  population  is  generated  randomly 
in  a  certain  range;  this  contains  many  individuals  and  is  used  to 
calculate  the  objective  functions  of  each  optimization  problem. 
Secondly,  the  individuals  are  ranked  using  a  fast  nondominated 
sorting  approach  according  to  their  objective  functions,  thus  the 
ranks  and  crowding-distances  are  obtained.  Thirdly,  a  binary 
tournament  selection  operator  is  used  to  choose  some  better  in¬ 
dividuals  as  parents.  Fourthly,  the  parents  are  crossed-over  and 
mutated  to  create  an  intermediate  generation.  Finally,  the  new 
generation  is  ranked  and  the  process  repeats.  The  Pareto  front  is 
obtained  by  selecting  those  individuals  whose  rank  value  is  equal  to 
one.  For  further  background  on  NSGA-II,  the  reader  is  referred  to 
Ref.  [24]. 

Fig.  4  shows  the  flowchart  of  the  modified  NSGA-II.  Here  we 
make  a  slight  change  to  the  crossover  and  mutation  process.  In  the 
original  NSGA-II  code  [25]  these  two  processes  are  simple  and 
inefficient.  Here  we  use  the  real  mutation  operator  and  discrete 
recombination  operator  in  the  Genetic  Algorithm  Toolbox  for 
MATLAB  26]  instead  of  those  in  the  original  NSGA-II  code.  The 
intermediates  generated  by  these  two  operators  have  higher  di¬ 
versity,  and  the  search  speed  is  higher. 

The  modified  NSGA-II  has  very  high  universality.  It  transforms 
the  optimization  problems  into  the  objective  function,  so  that  the 
optimization  problem  is  independent  of  the  genetic  and  non¬ 
dominated  sorting  process.  This  algorithm  is  very  suitable  for 
solving  highly  complex  and  large-scale  problems  such  as  parameter 
identification  of  a  multi-physics  model  for  Lithium-ion  batteries. 


4.  Parameter  identification  using  multi-objective  genetic 
algorithm 


4.2.  Parameter  identification 


4 A.  A  modified  NSGA-II 

The  Nondominated  Sorting  Genetic  Algorithm  (NSGA-II)  is  a 
multi-objective  optimization  approach  first  proposed  by  K.  Deb 
[24].  It  provides  a  set  of  optimal  solutions  known  as  a  Pareto  front. 


Fig.  2.  Cross  section  of  a  cylindrical  Li-ion  battery  [20]. 


Our  previous  work  [27  has  clarified  the  critical  parameters  to  be 
identified  and  their  search  range.  Some  of  these  25  parameters 
(such  as  De,  Rfiim,  etc.)  may  change  their  values  during  the  aging 
process,  and  therefore  the  identification  result  can  be  used  for 
researching  the  SOH  of  a  battery.  Other  parameters  (such  as  /<s,  E, 
etc.)  must  be  identified  because  their  values  are  extremely  difficult 
to  measure.  In  addition,  two  temperature-related  parameters,  the 
average  thermal  conductivity  in  the  radial  direction  Ar,  and  the  heat 
transfer  coefficient  h  are  also  included  in  this  analysis.  The  search 
ranges  of  these  parameters  are  shown  in  Table  4,  considering  both 
the  benchmark  values  from  references  and  the  probable  variation 
trend  during  the  aging  process.  For  temperature-related  parame¬ 
ters,  the  reference  value  at  the  reference  temperature  (298  K)  and 
the  activation  energy  in  the  Arrhenius  law  are  both  identified.  Some 
Li+  concentration-related  parameters  (i.e.  De  and  /ce)  are  treated  as 
constants  here,  the  reasons  being: 

(1)  There  is  no  particular  formal  relationship  between  the  Li+ 
concentration  and  the  parameter  values;  if  some  universal 
dependence  (e.g.  polynomial,  exponential  or  power  law) 
were  used,  it  would  incorporate  many  unknown  parameters 
to  be  identified,  making  the  identification  work  more  com¬ 
plex  or  even  impossible. 

(2)  In  Refs.  [10,19  ,  particular  equations  involving  electrolyte 
parameters  De  and  Ke  are  presented.  However,  during  the 
aging  process,  the  properties  of  the  electrolyte  change 
significantly,  as  shown  in  Refs.  [3,5  ,  and  these  equations  no 
longer  represent  the  real  values  of  those  parameters. 
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*(s) 


Fig.  3.  Current  profiles  of  (a)  LiFeP04  and  (b)  LiCo02  battery  in  simplified  DST. 


(3)  From  the  model  simulation  we  found  that  the  Li+  concen¬ 
tration  of  the  electrolyte  is  relatively  constant  during  the 
discharging  process  and  has  little  effect  on  the  parameter 
values. 

As  a  result,  we  treated  the  concentration-related  parameters  as 
constants  in  this  paper.  In  our  further  work,  these  parameters  will 
be  identified  in  different  aging  stages  and  used  for  the  analysis  of 
electrolyte  degradation. 

In  addition,  some  geometric  parameters  (e.g.  the  thickness  of 
electrode  layer)  and  material  properties  (e.g.  density,  specific  ca¬ 
pacity)  are  required  to  complete  the  model  simulation.  These  pa¬ 
rameters  have  been  obtained  from  battery  specifications,  literature 
values  or  direct  measurement  after  disassembling  real  batteries. 

In  this  article,  the  optimal  solution  of  an  identification  problem 
is  a  set  of  parameters,  with  which  the  simulated  data  should  pro¬ 
vide  a  best  fit  to  the  experimental  data.  So,  the  objective  functions 


Fig.  4.  Flowchart  of  the  modified  NSGA-II. 


are  defined  as  the  Sum  of  Squares  Error  (SSE)  between  the  model 
outputs  and  the  experimental  data  from  the  simplified  DST  at  15  °C 
and  30  °C.  The  mathematical  formulations  are  as  follows: 


N 


ObjF}  —  'y  ^  f^app,i,15°C  b/app  j  15° c 


i=l 

N 


minimizes 


ObjF2  —  ^  (rsh,i,15  -c  “  ^sh,i,15  °C J 
i=  1 

*  /  ~  v 

Obj F3  =  ^  ( 

b^app,i,30  °C  ^app,i,30  °C J 


(19) 


i=l 

N 


Objf4  —  ^  ^Tsh  j  30  «c  ^sh,i,30°c) 
i=l 


where  Uapp,i  and  Tshj  are  the  model  output  terminal  voltage  and 
surface  temperature  data,  Uappi  and  Tshi  are  the  experimental 
data,  and  N  is  the  number  of  data  points. 

The  objective  functions  are  minimized  by  using  NSGA-II.  As 
these  functions  approach  zero,  better  agreement  is  obtained  be¬ 
tween  the  model  output  and  the  experimental  data.  The  identifi¬ 
cation  code  is  implemented  in  MATLAB.  We  set  the  population  size 
NIND  =  400,  and  the  maximum  generation  MAXGEN  =  200.  Since 
the  simulation  time  of  the  simplified  DST  is  about  12  s  using  the 
multi-physics  model,  and  the  model  will  be  called  about  80,000 
times  during  the  whole  identification,  it  is  computationally 
expensive  and  the  total  time  consumption  is  unacceptable.  In  order 
to  speed  up  the  identification,  we  built  a  computing  cluster  with  20 
cores  using  an  HP  ProLiant  DL160  G6  server  (2x  Intel  Xeon  L5520  @ 
2.26  GHz  CPU,  16  GB  RAM)  and  five  PCs  (Intel  Core  i3-530  @ 
2.93  GHz  CPU,  4  GB  RAM)  based  on  the  MALTAB  Distributed 
Computing  Toolbox.  The  GA  identification  code  was  parallelized 
and  successfully  accelerated  by  a  factor  of  over  14  times  in  our 
previous  work  [28].  The  same  method  is  used  in  NSGA-II  here,  so 
the  actual  identification  time  is  about  19  h  for  each  battery. 

When  the  parameter  identification  is  complete,  NSGA-II  pro¬ 
vides  a  set  of  optimal  solutions  called  Pareto  fronts  for  each  type  of 
battery,  instead  of  a  single  optimal  solution.  The  individuals  on  the 
Pareto  front  represent  different  parameter  sets,  each  of  them  hav¬ 
ing  four  objective  function  values,  a  rank  value  and  a  crowding- 
distance  value.  Theoretically,  an  ideal  solution  should  exist  for 
which  all  objectives  (ObjFj  ~  ObjF4)  are  minimized.  However,  by 


Table  3 

Attribute  matrix. 


Individual 

Objective  functions 

F\ 

f2 

h 

f4 

Pi 

T\  1 

*21 

P31 

U 1 

P2 

O2 

*22 

*32 

J42 

Pj 

rV 

r2j 

r3 j 

r4j 

Pm 

r-im 

r2m 

r2m 

r4m 
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Table  4 

Search  range  and  the  final  identification  results  for  two  types  of  Li-ion  battery. 


Parameter 

Unit 

LiFeP04 

LiCo02 

Lower 

Upper 

Result 

Lower 

Upper 

Result 

Rs,a 

(pm) 

1 

9 

1.0247 

5 

15 

7.8502 

Rs,c 

(pm) 

0.01 

0.1 

0.0331 

5 

15 

11.2368 

fi^s.a.ref 

(m2  s’1) 

5  x  10'15 

1  x  10'13 

9.8910  x  10'14 

1  x  10'14 

9  x  10'14 

6.8814  x  10'14 

fi^s.c.ref 

(m2  s'1) 

5  x  10'19 

1  x  10  17 

4.2958  x  10  18 

5  x  10'15 

5  x  10'14 

1.8062  x  10'14 

^s,c 

(S  m-1) 

0.1 

5 

0.4977 

1 

20 

1.1901 

*0 

(-) 

0.75 

0.85 

0.7514 

0.65 

0.78 

0.6515 

yo 

(-) 

0.02 

0.08 

0.0367 

0.4 

0.5 

0.4347 

^s,a 

(-) 

0.5 

0.6 

0.5438 

0.3 

0.5 

0.4319 

£s,C 

(-) 

0.4 

0.6 

0.4135 

0.4 

0.65 

0.5610 

^film.a 

(Q  m2) 

0.002 

0.01 

0.0045 

0.001 

0.005 

0.0002 

Ce 

(mol  m'3) 

800 

1200 

1073.6 

900 

1200 

959.6 

fi^e.ref 

(m2  s'1) 

1  x  10'11 

1  x  10  9 

3.9378  x  10  11 

2  x  10'10 

1  x  10'9 

2.0241  x  10  10 

^e,ref 

(S  nr1) 

0.1 

2.5 

0.1917 

0.2 

4 

0.2101 

Rext 

(Q  m2) 

0.002 

0.005 

0.0031 

0.001 

0.005 

0.0038 

£e,a 

(-) 

0.3 

0.4 

0.3442 

0.3 

0.5 

0.3362 

£e,s 

(-) 

0.4 

0.6 

0.4068 

0.5 

0.6 

0.5081 

£e,c 

(-) 

0.3 

0.4 

0.3015 

0.3 

0.5 

0.3117 

^s.a.ref 

(m2-5  mol-0-5  s'1) 

CN 

1 

o 

y-* 

X 

y—* 

1  x  10  10 

3.9530  x  10  11 

1  x  10'11 

9  x  10'11 

6.9020  x  10'11 

^s.c.ref 

(m2-5  mol'0  5  s'1) 

1  x  10'13 

1  x  10  11 

7.8816  x  10  12 

1  x  10'12 

9  x  10'12 

7.8980  x  10  12 

H<s,a 

(KJ  mol-1) 

10 

50 

26.81 

10 

50 

39.57 

ERS,c 

(KJ  mol'1) 

10 

50 

31.92 

10 

50 

37.48 

EDs,a 

(KJ  mol'1) 

10 

50 

10.37 

10 

50 

42.77 

EDs,c 

(KJ  mol-1) 

10 

50 

12.86 

10 

50 

18.55 

EDe 

(KJ  mol'1) 

10 

50 

24.94 

5 

40 

37.04 

Ene 

(KJmor1) 

10 

50 

41.88 

10 

50 

34.70 

A 

(W  ITT1  K'1) 

0.1 

2 

0.1517 

0.1 

3 

1.9386 

H 

(W  m  2  K'1) 

5 

50 

28.41 

5 

50 

21.45 

analyzing  the  final  Pareto  front,  we  found  that  such  an  ideal  solu¬ 
tion  does  not  in  fact  exist.  Solutions  on  the  Pareto  front  could 
typically  only  achieve  minimization  of  one  or  two  objectives  and  no 
one  solution  could  optimally  satisfy  all  objectives  simultaneously. 
This  means  that  there  is  no  ideal  parameter  set  with  which  the 
simulated  data  of  Dapp  and  Ts h  could  perfectly  fit  the  experimental 
data  of  the  simplified  DST  at  both  15  °C  and  30  °C. 

However,  an  optimum  solution,  which  is  closest  to  the  ideal 
solution,  should  exist  on  the  Pareto  front.  Therefore,  a  reasonable 
and  effective  method  is  required  to  select  it  from  the  Pareto  front, 
and  this  chosen  solution  will  be  treated  as  the  final  identification 
result. 

4.3.  Selection  of  the  final  identification  result 

The  purpose  of  this  section  is  to  describe  the  procedure  for  se¬ 
lection  of  a  final  identification  result  from  each  Pareto  front  using 
the  values  of  objective  functions.  The  objectives  cannot  be 
compared  directly  because  of  their  different  magnitudes  and  units. 
This  is  a  typical  Multiple  Criteria  Decision  Making  (MCDM)  prob¬ 
lem,  for  which  the  Technique  for  Order  Preference  by  Similarity  to 
Ideal  Solution  (TOPSIS)  [29]  can  be  used  to  solve  it.  According  to  Eq. 
(19),  the  ideal  solution  of  the  parameter  identification  problem  is 
the  one  that  can  render  all  four  objective  functions  equal  to  zero,  or 
failing  that  equal  to  minimum  values. 

In  this  work,  the  four  values  of  the  objective  functions  of  each 
individual  compose  an  attribute  matrix.  For  the  LiFePCH  and  LiCo02 
batteries,  the  Pareto  fronts  each  have  400  solutions.  In  other  words, 
all  individuals  in  the  final  population  are  on  the  Pareto  front,  so  the 
dimension  of  the  attribute  matrix  is  (400, 4).  As  shown  in  Table  3,  pj 
( j  =  1,  2,  ...,  m)  stands  for  the  jth  individual  on  the  Pareto  front, 
Fi  =  (rn,  r,2,  ...,  V[ m)T  (i  =  1,2, 3, 4)  is  the  objective  function  vector  of 
the  ith  objective,  is  the  ith  objective  function  value  of  the  jth 
individual,  m  is  the  number  of  individuals  in  the  attribute  matrix. 


Since  the  rank  values  of  the  individuals  in  the  Pareto  front  are  all 
the  same  and  equal  to  1,  this  metric  is  excluded  from  further  study. 

Firstly,  the  attribute  matrix  is  reprocessed  according  to  the 
crowding-distance.  Crowding-distance  is  used  to  describe  the  dis¬ 
tance  between  different  individuals  in  NSGA-II.  When  the 
crowding-distance  of  an  individual  is  equals  infinity,  it  is  located  on 
the  edge  of  the  Pareto  front,  and  in  this  case  only  one  objective  is 
minimized.  Therefore,  the  individuals  whose  crowding-distance 
are  equal  to  infinity  are  removed  from  the  attribute  matrix. 

Secondly,  the  decision  matrix  is  normalized.  It  is  obvious  that 
each  element  r- •  in  the  decision  matrix  is  in  the  range  [0,1],  and  in 
the  vector  F/  the  normalized  attribute  of  the  optimal  individual 
(which  has  the  minimum  objective  function  value)  is  1,  while  the 
normalized  attribute  of  the  worst  individual  (which  has  the 
maximum  objective  function  value)  is  0.  The  normalization 
method  is: 


(20) 


where  F^*  and  F^  are  the  maximum  and  minimum  values 
respectively  in  vector  F/,  and  Ffpfi  =  ry  here. 

Thirdly,  we  define  an  ideal  solution  Z*  =  (1,1, 1,1 ),  and  a  negative 
ideal  solution  Z~  =  (0,0, 0,0).  Then  the  distance  between  an  indi¬ 
vidual  and  the  ideal  solution  is  obtained  as  follows,  using  the 
Euclidean  norm: 


E  ( 4  - 1)2 

i=  1 

Similarly,  the  distance  between  an  individual  and  the  negative 
ideal  solution  is: 


j=  1,2, ...,m 


(21) 
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Fig.  5.  Experimental  data  of  (a)  Ua pp,  (b)  TSh  for  LiFeP04  and  (c)  Uapp,  (d)  7Sh  for  LiCo02  battery  in  simplified  DST. 


Fig.  6.  Comparison  of  simulated  data  with  experimental  data  with  four  identification  objectives  for  the  LiFeP04  battery. 


E  (zv  ~  °)2 

i=l 

where  Zy  is  the  weighted  element  in  the  decision  matrix: 


j  =  1,2, ...,m 


(22) 


(i  =  1,2,3,4j  = 


1,2, 


(23) 


Sr  +S 

]  j 


j=  1,2 


(24) 


where  0  <q<  1.  The  closer  the  C*  is  to  1,  the  closer  an  individual 
is  to  the  ideal  solution. 


where  W[  is  the  weight  coefficient  of  the  ith  objective.  In  this  work 
W  =  (1,1, 1,1),  which  means  that  the  four  identification  objectives 
have  equal  importance  in  the  final  decision. 

Fourthly,  the  relative  adjacent  degree,  which  is  ascertained  via 
the  distance  between  individuals  and  the  ideal  solution  and  the 
distance  between  individuals  and  the  negative  ideal  solution,  is 
obtained  as  follows: 


Table  5 

Percentile  errors  and  average  errors  for  LiFeP04  battery  parameter  identification. 


Objectives 

Errorioo% 

Error75% 

Error50% 

Error25% 

ErroraVg 

(a)  Uapp  at  15  °C 

154  mV 

16.5  mV 

12.9  mV 

8.7  mV 

14.8  mV 

(b)  Tsh  at  15  °C 

0.38  °C 

0.12  °C 

0.07  °C 

0.04  °C 

0.10  °C 

(c)  t/app  at  30  °C 

33.6  mV 

10.5  mV 

8.1  mV 

4.9  mV 

9.1  mV 

(d)  TSh  at  30  °C 

0.34  °C 

o 

CD 

O 

n 

0.06  °C 

0.02  °C 

0.10  °C 
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Fig.  7.  Discharging  characteristics  at  0.5C,  1C  and  2C  for  the  LiFeP04  battery  at  15  °C. 


Finally,  all  individuals  in  the  attribute  matrix  are  sorted  ac¬ 
cording  to  the  relative  adjacent  degree.  The  one  that  has  the 
maximum  C*  is  the  final  result  of  parameter  identification. 

5.  Results  and  discussion 

The  variation  of  potential  and  surface  temperature  with  time 
during  the  DST  is  shown  for  each  type  of  battery  in  Fig.  5.  The 
overpotential  of  both  batteries  at  15  °C  is  greater  than  that  at  30  °C, 
implying  that  during  the  discharge  process  the  terminal  voltage  at 
15  °C  is  lower,  and  vice  versa  during  the  charging  process.  This  is 
consistent  with  a  lower  diffusion  coefficient,  lower  electrochemical 
reaction  rate  and  higher  internal  resistance  at  the  lower  tempera¬ 
ture.  These  temperature-related  phenomena  can  be  accounted  for 
by  the  Arrhenius  law  in  the  multi-physics  model  (Eq.  (15').  In 
addition,  higher  internal  resistance  also  leads  to  higher  heat  gen¬ 
eration  rate,  which  causes  the  surface  temperature  at  15  °C  to  rise 
more  than  that  at  30  °C  during  the  discharge  process.  Traditional 
single-objective  identification  uses  the  Dapp  curve  at  a  certain 
ambient  temperature  as  the  only  optimization  objective,  and 
therefore  the  identification  results  cannot  reliably  predict 
temperature-related  phenomena  at  other  ambient  temperatures. 
Compared  to  single-objective  identification  using  only  one  test,  the 
multi-objective  method  can  provide  more  accurate  and  robust 
results. 

Final  identification  results  for  the  two  batteries  are  listed  in 
Table  4.  Significantly,  for  the  26650-type  LiFePC^  battery,  the 
identified  value  of  Ar  is  0.1517  W  m_1  K-1,  which  is  very  close  to  that 
(0.15  ±  0.01  W  m-1  K-1)  measured  with  a  novel  experimental 
technique  in  Ref.  [21  .  This  demonstrates  that  the  thermal  imped¬ 
ance  model  presented  in  Section  2.3  is  accurate  and  effective,  and 
that  the  battery  surface  temperature  is  very  suitable  for  identifying 
thermal  parameters.  Furthermore,  the  identified  value  of  h  is 
20-30  W  nrr2  K-1  for  both  batteries,  much  higher  than  that  of  free 
convection  air  cooling.  Forced  air  convection  actually  occurs  on  the 
battery  surface,  because  the  thermal  chamber  uses  flowing  air  to 
control  the  ambient  temperature. 


The  accuracy  of  the  final  identification  results  will  be  discussed 
separately  for  each  type  of  battery  in  the  next  two  sections.  Vali¬ 
dation  studies  are  also  presented  by  comparing  the  simulated  and 
experimental  data  under  several  operating  conditions. 

5.1.  LiFeP04  battery 

Fig.  6  compares  the  simulated  data  with  the  experimental  data 
with  four  identification  objectives  for  the  LiFeP04  battery.  All 
simulated  curves  agree  well  with  the  experimental  curves,  the 
percentile  errors  and  average  errors  are  shown  in  Table  5.  Fig.  6(a) 
and  (b)  shows  the  Dapp  and  Ts h  comparison  curves  at  15  °C,  and 
Fig.  6(c)  and  (d)  shows  the  Dapp  and  Ts h  comparison  curves  at  30  °C, 
respectively. 

Discharging  tests  at  several  different  rates  at  15  °C  and  30  °C 
were  used  to  validate  the  accuracy  and  robustness  of  the  final 
identification  result.  The  terminal  voltage  and  surface  temperature 
data  are  compared  with  the  simulated  data  in  Figs.  7  and  8 
respectively.  The  red  (in  the  web  version)  marker  “O”  is  experi¬ 
mental  data,  and  the  blue  solid  line  is  the  model  simulation  using 
the  final  identification  result.  The  percentile  errors  and  the  average 
errors  between  the  simulated  and  experimental  data  are  also 
shown  in  each  sub-figure. 

It  is  observed  that,  under  all  operating  conditions,  the  simulated 
Dapp  curves  fit  the  experimental  data  very  well,  and  the  average 
errors  in  Dapp  never  exceed  21  mV  for  any  of  the  validation  condi¬ 
tions.  Similarly,  the  average  errors  in  Ts h  under  all  conditions  are 
below  0.19  °C.  Considering  the  temperature  measurement  resolu¬ 
tion  of  our  equipment  is  only  0.1  °C,  the  simulated  Ts h  curves  are 
also  in  good  agreement  with  the  experimental  data.  This  provides 
confidence  in  the  final  identification  result  and  the  robustness  of 
the  multi-physics  model  in  predicting  electrochemical  and  thermal 
behavior  of  real  batteries  under  different  operating  conditions. 

At  low  rate  discharge  process,  the  surface  temperature  curves 
present  as  wave-shape,  and  the  simulated  curves  seem  not  well  fit 
to  experimental  data.  Before  70  min  at  0.5C  discharging,  the  en¬ 
tropy  coefficient  of  the  whole  battery  changes  around  zero,  leading 
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Fig.  8.  Discharging  characteristics  at  0.5C,  1C  and  2C  for  the  LiFeP04  battery  at  30  °C. 
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Fig.  9.  Comparison  of  simulated  data  with  experimental  data  with  four  identification  objectives  for  the  LiCo02  battery. 


a  wave  of  the  surface  temperature.  And  the  total  temperature  rise  is 
only  1.1  °C  in  Fig.  7  and  0.7  °C  in  Fig.  8  for  0.5C  discharging, 
compared  to  the  equipment  resolution  (0.1  °C),  the  experimental 
data  is  actually  very  rough.  But  considering  the  average  error  is 
almost  equal  to  the  resolution,  the  relevance  between  the  simu¬ 
lated  curves  and  the  experimental  data  is  acceptable.  In  high  rate 
discharge  process,  the  ohmic  heat  which  is  almost  constant  be¬ 
comes  the  main  heat  generation  source  and  determines  the  bat¬ 
tery's  thermal  behavior,  and  the  total  temperature  rise  is  more  than 
2  °C,  so  the  Tsh  curves  look  more  smooth. 

5.2.  UC0O2  battery 

Similarly,  Fig.  9  compares  the  simulated  data  with  the  experi¬ 
mental  data  with  four  identification  objectives  for  the  LiCo02  bat¬ 
tery.  Fig.  9(a)  and  (b)  shows  the  Dapp  and  Ts h  comparison  curves  at 


15  °C,  and  Fig.  9(c)  and  (d)  shows  the  Dapp  and  Ts h  comparison 
curves  at  30  °C,  respectively.  The  percentile  errors  and  average 
errors  are  given  in  Table  6.  Again,  good  agreement  is  observed 
between  simulated  and  experimental  data. 

Discharging  tests  at  three  different  rates  (0.5C,  1C  and  1.5C)  at 
15  °C  and  30  °C  were  also  used  to  validate  the  final  identification 


Table  6 

Percentile  errors  and  average  errors  for  LiCo02  battery  parameter  identification. 


Objectives 

Errorioo% 

Error75% 

Error50% 

Error25% 

Erroravg 

(a)  Uapp  at  15  °C 

136  mV 

26.4  mV 

12.6  mV 

6.2  mV 

17.2  mV 

(b)  Tsh  at  15  °C 

0.78  °C 

0.28  °C 

0.21  °C 

0.11  °C 

0.23  °C 

(c)  Uapp  at  30  °C 

76.3  mV 

14.3  mV 

6.6  mV 

3.7  mV 

10.0  mV 

(d)  Tsh  at  30  °C 

0 

bo 

O 

n 

0.39  °C 

0.20  °C 

0.11  °C 

0.26  °C 
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Fig.  10.  Discharging  characteristics  at  0.5C,  1C  and  1.5C  for  the  LiCo02  battery  at  15  °C. 


results.  Comparisons  of  simulated  and  experimental  curves  are 
shown  in  Figs.  10  and  11.  At  30  °C,  the  simulated  Happ  curves  fit  the 
experimental  data  well,  and  the  maximum  error  never  exceeds 
64  mV.  But  at  15  °C,  there  are  some  discrepancies  between  the 
simulated  Happ  curves  and  the  experimental  data  during  the  latter 
part  of  the  discharge  process,  leading  to  a  maximum  error  of  over 
300  mV  at  high  discharge  rates.  Since  the  electrochemical  and 
thermal  parts  of  the  model  are  coupled,  an  inaccurate  voltage 
output  affects  the  temperature  output  and  vice  versa.  Both  the 
simulated  Happ  and  Ts h  curves  deviate  from  the  experimental  data 
in  that  period. 

The  discrepancy  is  more  apparent  at  higher  discharge  rates 
and  lower  ambient  temperature;  we  speculate  that  the  possible 
reasons  may  be:  (1)  The  properties  of  the  active  material  might 
be  changed  under  these  conditions,  and  the  related  parameters 
no  longer  constant  as  assumed  in  the  model.  (2)  The  model  ig¬ 
nores  some  phenomena,  such  as  side  reactions.  Model  simulation 
cannot  predict  real  battery  behavior  in  these  extreme  conditions. 
(3)  Some  design  and  geometrical  parameters  used  in  the  model 
may  be  inaccurate.  Such  parameters  are  determined  by  mea¬ 
surement  or  estimation  after  disassembling  a  real  LiCo02  battery. 
By  contrast,  for  the  LiFePC^  battery,  these  parameters  are  ob¬ 
tained  from  the  literature  [18,19],  which  might  imply  that  they 
are  more  accurate  and  have  been  experimentally  validated  by 
several  researchers.  (4)  The  model  entropy  coefficients  and 
reference  equilibrium  potential  might  be  different  from  that  in 
the  real  battery. 

However,  the  simulated  Happ  and  Ts h  curves  with  the  final 
identification  result  capture  the  major  characteristics  for  all  vali¬ 
dation  data,  such  as  the  plateau  of  terminal  voltage  and  the  trend  of 
surface  temperature.  Though  some  discrepancies  exist  in  compar¬ 
ison  to  the  LiFePCH  battery,  the  correctness  and  robustness  of  the 
final  result  for  the  UC0O2  battery  have  been  validated. 

By  the  way,  the  surface  temperature  evolution  of  these  two 
types  of  batteries  shows  varying  behavior  under  the  same  oper¬ 
ating  conditions,  because  the  entropy  coefficients  of  the  LiCo02  and 


LiFePCH  electrodes  are  very  different.  This  phenomenon  can  be 
used  to  distinguish  the  electrode  material  of  a  Lithium-ion  battery 
in  other  applications. 

In  order  to  increase  the  accuracy  of  parameter  identification  of 
the  multi-physics  model  using  NSGA-II  and  TOPSIS  method,  further 
work  will  be  carried  out  on  the  following  aspects: 

(1)  Improvements  to  the  accuracy  of  battery  surface  tempera¬ 
ture  would  increase  the  value  of  surface  temperature  as  an 
identification  objective. 

(2)  More  accurate  values  of  design  or  geometrical  parameters 
should  be  used  in  model  simulation,  in  order  to  increase  the 
identifiability  of  unknown  parameters.  If  such  accurate 
values  cannot  be  obtained  directly,  they  should  be  identified 
as  unknown  parameters. 

(3)  More  complex  operating  conditions  (i.e.  charge/discharge 
rates  and  ambient  temperatures)  should  be  used  as  identi¬ 
fication  tests  to  increase  the  robustness  of  the  identification 
result. 

(4)  In  the  TOPSIS  process,  a  reasonable  method  for  determining 
the  weighting  coefficients  should  be  developed  based  on  the 
credibility  of  identification  objectives.  For  example,  if  the 
accuracy  of  measured  Ts h  data  is  relatively  low,  the  objective 
weight  of  Tsh  in  Eq.  (23)  should  be  lower  than  that  of  other 
credible  objectives. 

6.  Conclusions 

In  this  work,  a  multi-objective  parameter  identification  method 
for  Lithium-ion  batteries  using  modified  NSGA-II  and  TOPSIS  al¬ 
gorithms  is  proposed.  Terminal  voltage  and  surface  temperature 
curves  at  15  °C  and  30  °C  are  used  as  four  identification  objectives. 
Two  Pareto  fronts  are  obtained  using  the  NSGA-II  algorithm  for  two 
types  of  real  batteries,  and  a  multiple  criteria  decision  making 
method  TOPSIS  is  used  to  select  the  final  identification  result  from 
the  Pareto  front.  Simulated  data  using  the  final  parameter 
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identification  results  are  in  good  agreement  with  the  experimental 
data  for  both  types  of  battery.  The  identified  parameters  can  indi¬ 
cate  the  battery's  internal  characteristics  accurately,  and  can  be 
used  for  analyzing  the  aging  mechanism  in  further  work.  It  should 
be  noted  that  the  model  outputs  with  the  final  results  for  the 
UC0O2  battery  show  some  discrepancies  from  the  experimental 
data  during  the  latter  part  of  the  discharge  process,  suggesting  that 
further  improvements  could  be  included  in  the  multi-objective 
parameter  identification  process. 
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Nomenclature 

Ce:  electrolyte  concentration  (mol  m-3) 

Cs:  concentration  of  Li+  in  the  intercalation  particle  (mol  m-3) 
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Ce/S:  concentration  of  Li+  at  the  surface  of  intercalation  particle  (mol  m-3) 
De :  electrolyte  diffusion  coefficient  (m2  s-1) 

Ds:  Li+  diffusion  coefficient  in  active  material  (m2  s-1) 

E:  activation  energy  (kj  mol-1) 

Eocv:  open  circuit  potential  of  electrode  (V) 

F;  Faraday's  constant  (=96487C  mol-1) 

h:  heat  transfer  coefficient  (W  m-2  K-1) 

iapp :  applied  current  (A) 

io:  exchange  current  density  (A  m-2) 

if.  solid  phase  current  density  (A  m-2) 

iy.  solution  phase  current  density  (A  m-2) 

is:  electrochemical  reaction  current  density  (A  m-2) 

ju:  wall  flux  of  Li+  on  the  intercalation  particle  of  electrode  (mol  m-2  s-1) 

ks:  electrochemical  reaction  rate  (m25  mol-0-5  s-1) 

q:  heat  exchange  rate  (W  m-2) 

Q:  heat  generation  rate  (W  m-3) 

R:  radius  of  the  cylindrical  battery  (m) 

Rext:  contact  resistance  (Q  m2) 

Rfiim:  SEI  film  resistance  (Q  m2) 

Rs:  intercalation  particle  radius  of  electrode  (m) 
tg  :  Li+  transference  number  in  the  electrolyte 
T:  battery  temperature  (K) 

Tam:  ambient  temperature  (K) 

TSh •  surface  temperature  of  the  battery  (°C) 


UapP:  applied  potential  (V) 
Greek  symbols 


a:  transfer  coefficient  of  the  electrochemical  reactions 
e:  emissivity 

es:  volume  fraction  of  active  material  in  electrode 
ip.  overpotential  (V) 

Ke:  ionic  conductivity  of  the  electrolyte  (S  m-1) 

Ar;  thermal  conductivity  in  the  radial  direction  (W  m-1 
o:  Stefan— Boltzmann  constant  (=5.6704  x  10-8  W  m-2 
as:  solid  phase  conductivity  (S  m-1) 

4>s:  solid  phase  potential  (V) 

4>e:  solution  phase  potential  (V) 

Subscript 


a:  anode 
c;  cathode 
s:  separator 
0:  initial  value 
ref:  reference  value 


